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Abstract. We consider the inverse conductivity problem in a strictly convex domain 
whose boundary is not known. Usually the numerical reconstruction from the measured 
current and voltage data is done assuming the domain has a known fixed geometry. How- 
ever, in practical applications the geometry of the domain is usually not known. This 
introduces an error, and effectively changes the problem into an anisotropic one. The 
main result of this paper is a uniqueness result characterizing the isotropic conductivities 
on convex domains in terms of measurements done on a different domain, which we call 
the model domain, up to an affine isometry. As data for the inverse problem, we assume 
the Robin-to-Neumann map and the contact impedance function on the boundary of the 
model domain to be given. Also, we present a minimization algorithm based on the use of 
Cotton-York tensor, that finds the pushforward of the isotropic conductivity to our model 
domain, and also finds the boundary of the original domain up to an affine isometry. This 
algorithm works also in dimensions higher than three, but then the Cotton- York tensor 
has to replaced with the Weyl-tensor. 

AMS classification: 35J25, 35R30, 58J32. 

Keywords: Inverse conductivity problem, electrical impedance tomography, unknown 
boundary, Cotton- York tensor. 

1. Introduction. We consider the electrical impedance tomography problem (EIT for 
short), i.e. the determination of the unknown isotropic conductivity distribution inside a 
domain in M 3 , for example the human thorax, from voltage and current measurements 
made on the boundary. Mathematically this is formulated as follows: Let Q be the mea- 
surement domain, and denote by 7 the bounded and strictly positive function describing 
the conductivity in Q. The voltage potential u satisfies in Q the equation 

(1.1) V-7Vw = 0. 

To uniquely fix the solution u it is enough to give its value on the boundary. Let this be /. 
In the idealized case, when the contact impedance of the measurement device is zero, one 
measures for all voltage distributions u\qm = f on the boundary the corresponding current 
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flux through the boundary, jdy/dv, where v is the exterior unit normal to dQ. Mathe- 
matically this amounts to the knowledge of the Dirichlet-Neumann map A corresponding 
to 7, i.e. the map taking the Dirichlet boundary values to the corresponding Neumann 
boundary values of the solution to (jl.ljl . 

du 

A : u qm i-> Itt- 
av 

The Calderon's inverse problem is then to to reconstruct 7 from A. The problem was orig- 
inally proposed by Calderon jH] in 1980 and then solved in dimensions three and higher for 
isotropic conductivities which are C°°-smooth in arid [21]. The smoothness require- 
ments have been since relaxed, and currently the best known result is |27] with unique 
determination of conductivities in W 3 ^ 2 '°°, see also JU] for a somewhat different approach 
to the lack of smoothness. In two dimensions the first global result is due to Nachman 
(|25|). and later Astala and Paivarinta showed in [3] that uniqueness holds also for general 
isotropic L°°-conductivities. For the corresponding anisotropic case, see |3J I2H1 I2H 122] and 
numerical implementations of the methods with simulated and real data, see [T31 UTT] 

Assuming that the measured Dirichlet-to-Neumann map A meas is given, an often used 
method to solve the EIT-problem is to minimize 

||A mca s - A CT || 2 + a|H|x 

for a defined in terms of some triangulation of Q and || • \\x is some regularization norm; 
here A a is the Dirichlet-Neumann map corresponding to the conductivity a. One then also 
fixes the geometry of Q by assuming that it is for example a ball or an ellipsoid. Now, if 
our measurements have no error a Bayesian interpretation of this problem as a search of 
an MAP-estimate suggests that a = 0. Usually, the given data A meas does not correspond 
to any isotropic conductivity in the model domain. The reason for this is that there is no 
conformal map deforming the original domain to the model domain. Therefore, in solving 
the minimization problem we obtain an incorrect solution a. This means that a systematic 
error in modeling causes a systematic error to the reconstruction. In particular, if we 
consider linearization 7 = 70 + £71 where 70 is given known background conductivity and e 
is small, it seems that a localized perturbation 71 gives a reconstruction a = 70 + e<J\ where 
the reconstructed perturbation <7i is not localized. This is clearly seen in brain-activity 
measurements, see [15] and [To] . 

This work is continuation of ^K] where the corresponding question in two dimensions was 
studied: we proved that on the model domain there is a unique (anisotropic) conductivity 
with minimal anisotropy. This follows from a result of Strebel saying that among all 
quasiconformal self-maps of the unit disk with a fixed boundary value there is a unique one 
with minimal complex dilation. In higher dimensions there are several new issues. First of 
all, the non-uniqueness due to anisotropy is not understood, except in the case when both 
the domain and the conductivity function are the real analytic ([221, EH])- Also, as we 
already mentioned, in the plane case one could use the theory of quasiconformal maps to 
break the non-uniqueness. The higher dimensional analogue of this is unknown. Finally, 
there is no analogue of the Riemann mapping theorem that we could use. 

The structure of this paper is the following. In the first part, consisting of sections 2-4, 
we present the uniqueness results we have on the problem. It is worth noting that we choose 
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to work with the Robin-to-Neumann (RN) map instead of the Dirichlet-to-Neumann (DN) 
map described above. Mathematically they are equivalent, as we will show, but the RN- 
map is a better model for the actual measurement configuration since it takes into account 
the contact impedances at dQ Also, we assume the the function modeling the contact 
impedances of the electrodes is known. There are two key ideas how we compensate for 
our lack of understanding of the full anisotropic problem. First thing is to note that if an 
isotropic conductivity is pushed forward by a diffeomorphism, the resulting conductivity is 
still conformally flat, and in three dimensions this is equivalent with the vanishing of the 
Cotton-York tensor. Secondly, we assume that our original domain is strictly convex, and 
then the Cohn-Vossen theorem can be used to determine the original boundary dfl up to 
rigid motions. 

In the second part we develop an algorithm for finding the shape of the domain Q and 
the conductivity inside using a minimization technique. Important feature is that we do 
not have to construct an embedding of the boundary to the Euclidean space. We plan to 
report on the numerical implementation of our algorithm in a separate article. 

2. Measurements. Let Q C M n , n > 3 be a strictly convex domain, and denote by 
7 = {Y"* ( x ))?j=i the symmetric real valued matrix describing the conductivity in Q. We 
assume that the matrix is bounded from above and from below, that is, for some C, c > 
we have 

(2.2) c||e|| 2 < (&7G«0O < C1KII 2 , for all xGSl. 

We will state the precise smoothness of 7 later. We start by consider the EIT problem 
with continuous boundary data. Instead of the Dirichlet-to-Neumann map we will use the 
Robin-to-Neumann map defined below that corresponds better to the measurements done 
in practice. We discuss later in this section the relation of the continuous model and the 
electrode measurements made in practice. 

For the electrical potential u we write the model 

(2.3) V-7Vw = 0, xett, 

(2.4) (zu- -fVu + u)\ dn = h, 

where h is the Robin-boundary value of the potential and z is a function describing the 
contact impedance on the boundary. The contact impedance models the impedance that 
is caused by electro-chemical phenomena at the interface of the skin and the measurement 
electrodes in practical measurements [7j. 

In mathematical terms, the perfect boundary measurements are modeled by the Robin- 
to-Neumann map R = R zn given by 

R : h 1— > v- 7Vw|an 

that maps the potential on the boundary to the current across the boundary. Next we 
relate this continuous model for measurements done in practice. 

The physically realistic measurements are usually modeled by the following complete 
electrode model (see |7J |33]): Let e,- C dQ, j = 1, . . . , J be disjoint open sets of the 
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boundary modeling the electrodes that are used for the measurements. Let u solve the 
equation 



(2.5) V-7Vw = infi, 

(2.6) ZjV-iSIv + v\ Cj = Vj, 
(2-7) "•7Vu|an\u/ =1 « y = °> 

where Vj are constants representing electric potentials on electrode Cj. Then, one measures 
the currents observed on the electrodes, given by 



1 f 

Ij = - — r / v- jVv(x) ds(x), j = 1, . . . , J. 

\ e j\ Jej 



Thus the electrode measurements are given by map E : M, J —>■ M. J , E(V\,...,Vj) = 
(ii, . . . , Ij). We say that E is the electrode measurement matrix for (dfl, 7, e±, . . . , ej, 
Zi, ■ ■ ■ , Zj). 

The complete electrode model can alternatively be defined as follows: the Robin-to- 
Neumann map R v is given by R v f = is-jVulgn where u is the solution of 

(2.8) V-7Vm = infi 

zv-^Wv + 7]v\ d n = h, 

where z G C°°(dQ) is such that its restriction to the electrode Cj is equal to the constant 
Zj and 77 = J2j=i Xe 3 , where x ej is the characteristic function of electrode tj. 

We associate to the electrode measurement matrix and to the complete electrode model 
also the corresponding quadratic forms E : M J xlR J — > M. and R v : H^^^dVl) xH^ 1 / 2 (d^l) — > 
K. given by 

(2.9) E[V 1 V] = Y](EV)jVj\e j l R v [h,h] = [ (R v h)hds. 

I 1 ' ' Jan 

These have the following simple relation to each other: Let S = span(x e ^ : j — 1, . . . , J) C 
H~ l / 2 {dn) and define M : V = {Vj) J j=l ^ ^ J =1 V jXej to be a map M : W -> S. Then 

(2.10) E[V, V] = R n [MV, MV\. 

By (12. the electrode measurement matrix can be viewed as the discretization of the 
form R v . By increasing the number of the electodes and making the gaps between them 
smaller, we can assume that rj — > 1. In this case R n approximates the Robin-to- Neumann 
map R ljZ . Note that E(V, V) corresponds to the power needed to maintain the voltages V 
in electrodes. 

In practical EIT experiments, one places a set of measurement electrodes on the bound- 
ary dQ, e.g., around the chest of the patient. All the traditional approaches to the numerical 
EIT reconstruction assume that the shape of the domain Q is known and the only unknown 
is the conductivity 7. However, in most EIT experiments the boundary of the body Q is 
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not known accurately and since there are no practically reliable measurement methods 
available for the determination of the boundary, the EIT image reconstruction problem is 
typically solved using an approximate model domain Q, which represents our best guess for 
the shape of the true body Q. However, it has been noticed that the use of slightly incor- 
rect model for the body O in the numerical reconstruction can lead to serious artefacts in 
reconstructed images fUJ • This situation is our paradigm for the EIT problem when 
the boundary is unknown. Next we analyze how the deformation of the domain affects 
measurements. 



3. Deformations of the domain. In this section we analyze the behavior of the elec- 
trode models under a diffeomorphism. Let's consider first the Robin-to-Neumann map R. 
The corresponding quadratic form, which we still denote by R, is given on the diagonal by 

(3.11) R[h,h]= / (u + zv- 7Vw)z/- 7VW cLSe = \ r )Vu-Vudx + I z\v- ^Vu\ 2 cLSe 

Jan Jn Jan 

where h E H^^^dVt), u solves (|2.8|h and (ISe is the Euclidean volume form (or area) 
of dVL. The value R[h, h] corresponds to the power needed to maintain the current h on 
the boundary. From the mathematical viewpoint, using the (incorrect) model domain Q 
instead of the original domain Q can be viewed as a deformation of the original domain. 
Thus, let us next consider what happens to the conductivity equation when the domain Q 
is deformed to Q. Assume that F : Q —>■ Q is a sufficiently smooth orientation preserving 
map with sufficiently smooth inverse F^ 1 : Q — > Q. Let / : dQ — > dQ be the restriction of 
F on the boundary. When u is a solution of V- 7VU = in Q, then u(x) = u(F~ l (x)) and 
h(x) = /i(/ _1 (x)) satisfy the conductivity equation 

(3.12) V-7W = 0, hi ft, 

ZV- 7VS + u\ 9 q = h. 

Here h(x) = h{f^ 1 {x)) 1 v is the unit normal vector of <9ft, z is the deformed contact 
impedance and 7 is the conductivity 



(3.i3) ~ l{x) =ny)ity)tny)Y 



detF'(y)| 



y=F~ 1 {x) 



where F' = DF is the Jacobian of the map F. Note that even if 7 is isotropic, i.e., scalar 
valued the deformed conductivity 7 can be anisotropic i.e., matrix valued. 

To determine the deformed contact impedance z, we consider the corresponding invariant 
(n — l)-form 

j := v^vudSE e n 71 - 1 (an) 

corresponding to the current flux through the boundary. Next we denote x = F(x). A 
straightforward application of chain rule gives that 

V-^Vu\ dh = ((detDF)-V Vu) of 1 !^ 
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since F was orientation preserving and DF is the Jacobian of F in boundary normal 
coordinates associated to the surface dQ C M n . In these coordinates det DF\q^ = detDf, 
where det Df is the determinant of the differential of the the boundary map / : dQ — > dQ. 
We note that (detDf o Z^ 1 )/*^^) = cISe, where cISe and cISe are Euclidean volume 
forms of dQ and dQ, respectively. Hence, zv ■ Vw transforms as an invariantly defined 
function when the contact impedance is interpreted as a density, i.e. 

(3.14) z(x) = (det Df(x))z(x) 

where f(x) = x. Now we see that the boundary measurements are invariant: When 
/ : dfl — > dQ is the restriction of F : Q — > Q, we say that the map R = f*R z<1 , defined by 

((UR zrt )h)(x) = (R Z)1 (hof))( y )\ y=f _ 1{x) , heH^ 2 (dn) 

has that property that R = Rs,^- We call R the push forward of R za by /. 

Is is also worth noting that in formula (|3.11j) the integral over Q, as well as the integral 
over the boundary are invariant because of the deformation rule (I3.14j) for the contact 
impedance z, that is, we have 

R[h 1 h'}=R[hof-\h'of-\ 

for h, ti G H-^ 2 (dQ). 

4. Uniqueness results Now we are ready to give the exact set-up of the problem 
we consider: We want to recover an image of the unknown conductivity 7 in Q from the 
measurements of Robin-to-Neumann map, and we assume a priori that 7 is isotropic. We 
assume z, dQ and R are not known. Instead, let Q, called the model domain, be our 
best guess for the domain and let f m : dfl — > dVL m be a diffeomorphism modeling the 
approximate knowledge of the boundary. 

As the data for the inverse problem, we assume that we are given the boundary of 
the model domain dQ, the function z o f~ l corresponding to the contact impedance of 
electrodes, and the Robin-to-Neumann map R = (f m )*R- Note the discrete analog of this 
data is to know the voltage-to-power form V 1— > E(V, V) and the contact impedances of 
the electrodes, but not the location of the electrodes or the boundary of the domain. It is 
reasonable to assume that the contact impedance z o f^ 1 on the boundary of the model 
domain is known since we can observe and set up the contact impedances of the electrodes 
the way we want. Hence we have on the boundary of our model domain dQ a boundary 
map R that does not generally correspond to any isotropic conductivity. Furthermore, 
we saw above that there are many anisotropic conductivities for which Robin-to-Neumann 
map is the given map R. Next we show that the existence of the "underlying" isotropic 
conductivity in Q gives the uniqueness in Q up to a diffeomorphism and that the domain 
Q and the isotropic conductivity on it can be uniquely determined. 

Theorem 4.1. Let Q C M n , n > 3 be a bounded, strictly convex, C°° -domain. Assume 
that 7 G C°°(Q) is an isotropic conductivity, z G C°°(dQ), z > be a contact impedance, 
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and i? 7jZ the corresponding Robin-to -Neumann map. Let fl be a model of the domain 
satisfying the same regularity assumptions as fl, and f m : dfl — > dfl be a C°° -smooth 
orientation preserving diffeomorphism. 

Assume that we are given dfl, the values of the contact impedance z(f~ 1 (x)), x G dfl 
and the map R = (f m )#Rj )Z . Then we can determine fl up to a rigid motion T and the 
conductivity 7 o T _1 on the reconstructed domain T(fi). 

We recall also that rigid motion is an affine isometry T : R" — > M. n . 

Proof. Assume we are given R and the values of the contact impedance, that is, 
the function Let F m : fl — > fl be an orientation preserving diffeomorphism 

satisfying F m \QQ = f m . As noted before, R = where z(x) = det(D f m ) z(f~ l (x)) is the 
contact impedance on dfl and 7 = (F m )*7 is the push forward of 7 in F m . The Robin-to- 
Neumann map is a classical pseudodifferential operator of order zero, with principal symbol 
1/z, and hence R determines z. Since we assume also zo / _1 known, we can determine the 
determinant det Df m ; note that this gives the change of boundary area under deformation 
f m . For the rest of the proof denote (3 = det Df m . Also, this implies that we can find the 
Dirichlet-to- Neumann map A^ = on dQ, that is, the map taking the Dirichlet 

boundary values to Neumann boundary values. 

The Riemannian metric corresponding to the isotropic conductivity 7 = j(x)I in Q is 
given by 

9ij (x) = det( 7 (x)/) 1 /( 2 -") (l(x)I)- 1 = 7(x) 2 ^-%. 

Then, if A g is the Laplace-Beltrami-operator corresponding to the metric g, we have 
A g = I £7 1 1 / 2 ^7 • 7V111 where \g\ = det(^y). This metric is an invariant object and in the 
deformation F m it is transformed to the metric g = (F m )*g in Q. By [22], the Dirichlet- 
to-Neumann map determines the metric tensor "g^ on the boundary in the boundary 
normal coordinates. In particular, if we consider dfl as a submanifold of lR n with the 
metric h = i*(g) inherited from (fl,g) where % : dfl —>■ fl is the identity map, we see 
that our boundary data determines the metric h on dfl. Let now metric h = i*(g) be the 
corresponding metric on dfl, where i : dfl —>■ fl is the identical embedding. Then we have 

(4.15) h=(f m ) m h, h = 1 2 ^h E 

where h E is the Euclidean metric of dfl. Denote by h E = (f m )*h E the metric tensor on 
dfl, i.e. the push-forward of the Euclidean metric of dfl by f m . Recall that dSs and dSs 
are the Euclidean volume forms of dfl and dfl, respectively. Then the Riemannian volume 
forms dS^ and dSh of the metrics h and h respectively satisfy 

dS~ h = (fmUdS h ) = lif-'imfmWSE) = ( 7 /3) o f~\x) dS E 

on dfl. As (3 was already determined, this shows that we can find l{f^ l l {x)), x G dfl and 
hence by (I4.15j) we can determine the metric 

h E = 1 (f-\x))-^h. 
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In other words, if we consider dQ as an abstract manifold that can be embedded to dQ C R n , 
we have found the metric tensor on dfl corresponding to the Euclidean metric of dQ. 
By the Cohn-Vossen rigidity theorem, intrinsically isometric C 2 -smooth surfaces that are 
boundaries of a strictly convex body are congruent in a rigid motion. For uniqueness, see 
e.g. [3U1 Thm. V and VI] and also (H1II21- This means that the boundary data determines 
uniquely the map T o jf" 1 , where T is a rigid motion. Hence we can find the surface T(dQ) 
and on it the map T^A^ = T*A 7 . Using the uniqueness of of the isotropic inverse problem 
1211, we see that the boundary data determines 7 o T _1 . □ 

Note that the construction of the surface dQ C R 3 from the intrinsic metric h E is a more 
delicate issue, see |261 I29| . hence we take care to avoid it. 

5. A reconstruction algorithm and the use of conformal flatness. In this section 
we consider the case n = 3, even though the considerations could be generalized for n > 4 
by changing the Cotton- York tensor to Weyl tensor in our considerations (see Appendix). 
As noted before, an actual construction of the isometric embedding of an abstract manifold 
to Euclidean space is complicated and thus we try to avoid it. 

We want to find an anisotropic conductivity 77 such that R-g^ = R assuming that R = 
{fm)*Rz,' r where 7 is an isotropic conductivity. Clearly, when F m : Q — > Q is diffeomorphism 
satisfying F m \on = f m , the anisotropic conductivity (F m )*7 is a solution of the inverse 
problem, but it is not unique. However, we also know that (-F m )*7 has a conformally 
flat structure and this fact will help in solving the inverse problem as we will see. Note 
that in principle, one could start to solve the inverse problem by minimizing over all pairs 
(Q, a) of smooth domains Q C R n and all isotropic conductivities a in Q. However, the 
minimization over domains is complicated, and our objective is to find a reasonably simple 
minimization algorithm where we minimize over conductivities in the fixed model domain 
Q with an appropriately chosen cost function. 

Let i] = (F TO )*7 be a possibly anisotropic conductivity in Q such that 7 is isotropic. As 
already noted it defines a Riemannian metric g on Q, given by 

From now on we will use the Einstein summation convention and omit the summation 
symbols. As i^" 1 : Q — > Q can be considered as coordinates, we see that in proper 
coordinates the metric g is a scalar function times Euclidean metric, that is, g is conformally 
flat. This means that 

9ij (x) = e~ 2 ^%(x) 

where ~g~ij(x) is a metric with zero curvature tensor (i.e. flat metric) and a(x) G R. By 
[Sj (for original work, see jH]), the conformal flatness of the metric g in three dimensions 
is equivalent to the vanishing of the Cotton- York tensor C = C^ corresponding to g (see 
Appendix). Note that we can choose a = 2z^log7 and g = (F m )*(8ij). By [HI formulae 
(28.18) and (14.1)], a satisfies a differential equation (with n = 3) 

(5.16) Oij = — RiCij + — -i- — g i:j R - - gijg lm oia k , i, j = 1, . . . , n 

n — 2 2(n — l)(n — 2) 2 



where Ric^ is the Ricci curvature tensor of g, R is the scalar curvature or g 1 and 

<9cr , a 

° k = dx*' a ' lj = l(Tj ~ ai ° j ' W &i = 

where V e . is the covariant derivative with respect to metric g. Thus if g is given, (j5.16|) is a 
second order nonlinear differential equation for a. By p. 92], the equations ()5.16|) satisfy 
the sufficient integrability conditions to be locally solvable if and only if the Cotton- York 
tensor vanishes. Note that the existence of the isotropic conductivity 7 in Q gives a solution 
for these equations. 

Consider now the following algorithm: 

Data: Assume that we are given dfl m , R = (f m )*R ltZ and z o f' 1 on dfl m . 

Aim: We look for a metric g corresponding to the conductivity 7 and 'z such that on <9fl m 
R = R^^ and z = (f m )*z. 

Algorithm: 

1. Determine the two leading terms in the symbolic expansion of R. They determine 
a contact impedance z and a metric 'g on dQ such that if R = Ry^ then z = z and 

i*(g) =9- 

2. Form the ratio of reconstructed i.e. z, and measured contact impedances 

r(x) := Z ^S, X) \ xedQ. 
z{x) 

Note that then 

r(x)(f m )*(dS E ) = dS E 
since the contact impedances transformed as densities. 

3. Let dSg be the volume form of "g on dfl. Then 

dS d = (det g) l/2 dS E - 

Define 

7 = (det g) 1/2 f. 

With this choice 7 will satisfy 7(x) = 7(/m 1 (^)) f° r ^ ^ 

4. Define the boundary value a for the function a by 

= log (j(x)) , 5; G dfl. 

2 — n 
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5. Solve the minimization problem 



min F T (z, er, 7) + a H(z, a, 7) 
where H(z, 7) is a regularization functional, say 
H(z,a,j) = \\z\\ H 8 ( n) + 117 



{si) "T" 11 n\ H 8{si) n a liH8(n)) 



|2 



« > is a regularization parameter, and 
Ft(z,&, t) = ll-R ~~ ^. ? il£(j?-i/2(an)) 

+ ll <^l))" f( ^ )llL2( ^ 
+riirn 2 - 

^' H^llL 2 (n) 

+ Il cr laf2 ~~ a \\h^{dsi) 



+ E K - (-^2 mC - + 2(n-l)(n-2) ^ R - 2^"^) 



|2 



where r > 0, g is the metric tensor corresponding to 7, C = Cy is the Cotton- York 
tensor of g, and finally Ric and i? are the Ricci curvature and scalar curvature tensors 
respectively of g. 

Note that above value of the Cotton- York tensor at x G Q, Cij(x), the Ricci curvature 
tensors Rij(x), and the scalar curvature R(x) depend on the values of the conductivity r\ 
and its three first derivatives at x. 

Proposition 5.1. Let Q c M 3 be a bounded, strictly convex, C°° -domain. Assume that 
7 G C°°(Q) is an isotropic conductivity, z G C°°(dQ), z > be a contact impedance 
and R JjZ be the corresponding Robin-to-Neumann map. Let Q be a model of the domain 
satisfying the same regularity assumptions as Q, and f m : dQ — > dQ be a C°° -smooth 
diffeomorphism. 

Assume that we are given dQ, the values of the contact impedance x G dQ, 

and the map R = (f m )*R-y, z . 

Let t > 0. Then minima of F T (z,a,^f) is zero and any minimizers z, a and 7 of 
F T (z, a, 7) satisfy z = (f m )* z and there is a diffeomorphism F : Q — > Q such that F\g n = 
f m , 7 = F*y and a = - log 7. 

Proof. Assume first that r > 0. The minimizer exists because of existence of fi, 7,2 
and a, and the minimum is zero. Let z, a and g be some minimizers of F T . As then the 
Cotton- York tensor is zero and the equations f)5. 16j) are valid, if follows from jU], that the 
metric g^ = exp(2a(x))gij(x) , x G Q is flat. Since = R, we have z = (f m )*z and 
the metric g corresponding to 7 has to satisfy on the boundary i*g = g. This, and the 
vanishing of F T imply that 

i*g = exp(2a)i*g = exp(2a)g = exp(2a)(f m )*(jh E ) = 
= exp(2a)j(f m ) :¥ (h E ) = (f rn )*(h E ). 
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Consider now (fl, g) as a Riemannian manifold. As g is flat, we know that (Q,g) can be 
embedded isometrically to domain Q C M ra . Let : — > Q be this embedding. Since 
i*9 = {fm)*{h>E), h follows from the Cohn-Vossen rigidity theorem that the boundary dQo 
and dQ are congruent in a rigid motion T and k o f m = T\qq. Then (T _1 o &;)*7 is isotropic 
conductivity, the contact impedances of (T _1 o k)*z and z coincide, and the Robin-to- 
Neumann maps of (T _1 o k)*cr and a coincide. By the uniqueness of the isotropic inverse 
conductivity problem [33], (T _1 o k),/y = 7. This proves the claim in the case r > 0. 

Next, consider the case r = 0. Again, minimizer exists because of existence of Q,j,z 
and a, and the minimum is zero. Let z, a and g be some minimizers. As the minima of F T 
is zero, the equations (|5.16|) are valid. By [§1 p. 92], the solutions a satisfy the integrability 
conditions 

(5.17) VfcO-jj - VjCTik = criR l ijk , i,j,k — l,...,n 

that imply that the conformal covariant satisfies Rijk vanishes. Thus the Cotton- York 
tensor Cy is zero. This means that the minimizers z, a and g of F T with r = are also 
minimizers of F T with any r > 0. □ 

One can think of r as a regularization parameter: in general the solvability properties 
of equations ()5.16|) are sensitive to the compatibility conditions, i.e. the vanishing of the 
Cotton- York (or the Weyl tensor in higher dimensions). 

To find the domain Q, we can continue the above algorithm by applying the fact that 
conformally Euclidean manifold of dimeision n can be conformally embedded to W 1 in a 
constructive way (cf. [E]). 

6. In the previous steps 1.-5. we have found metric tensors g and g = e 2cr g on Q such 
that lj = F*(g) and g = F*(g E ) where g is the metric corresponding to the metric 7 
on f2, g E is the Euclidean metric on Q, and F : Q — > Q is some diffeomorphism. 

Let y 6 Q and find geodesies ~p y ^(s) starting from y with respect to the metric ~g. 
We parametrize these geodesies in such a way that JI y ^(0) = y and d s Jiy^(0) = £ is a 
unit tangent vector of the tangent space (T y Q,g). These geodesies correspond to the 
halflines in M. 3 starting from some point yo 6 fi. Let J : (T y Q,g) — > IR 3 be a linear 
isometry and define a map k : Q — > M. 3 by setting 

Then k o F : f2 — > M 3 is an afline isometry that extends to a rigid motion T : M 3 — > M. 3 
with = 0. Thus we can find k(Q) = T(Q), = TVy, and = T*^ 

Thus we have shown the followig reconstruction result. 

Corollary 5.2. Let n,j,z,Q and f m be as in Proposition I5.il Assume that we are 
given dQ, the contact impedance ^(/~ 1 (x)), x G dQ, and the Robin-to- Neumann map 
R — (fm)*R^,z- Then the algorithm 1.-6. determines Q, 7, and z up to a rigid motion 
T : M 3 -> M 3 . ' 

We intend to investigate the numerical implementation of the method and give numerical 
test results in part II of this paper. 



11 



Appendix. Here we define the conformal curvature tensors. We say that a metric in 
a domain Q C W 1 is conformally flat if there is a scalar function a(x) > such that the 
curvature of tensor of a(x)gij(x) is identically zero. 

First, let 7 be an isotropic conductivity, i.e., a smooth positive function in Q and F : 
Q — > Q be a diffeomorphism. Let 77 = F*7 be a possibly anisotropic conductivity in fl. It 
defines a Riemannian metric g on Q, given by 

As : f2 — ► can be considered as coordinates, we see that in proper coordinates the 
metric g is a scalar function times Euclidean metric, that is, g is conformally flat. 

Next we consider a general metric tensor g^ and recall facts concerning its conformal flat- 
ness Note that below we use the Einstein summation convention and omit the summation 
symbols when possible. The following tensors are related to conformal flatness. 

(a) Assume that n = 3. Then the conformal covariant, given in terms of curvature tensors 
(see the explanation on notations below) is 

Rijk = Vfc-Rjj — VjRik + _ i\ (dik^jR — gij^kR)- 

In the three dimensional case R^k defines a tensor that can be considered as vector 
valued 2- form Rij^dx^ A dx k . Operating with the Hodge operator * to this 2-form, 



we obtain the Cotton- York tensor 



Cij = g kp g lq V k (Ru-±Rgu)e pqj , 



where e pq j is the Levi-Civita permutation symbol, 
(b) Assume that n > 4. Then the Weyl tensor is 

WijM = RijM + — — -{guRkj + gjkRkj — gtkRij — gjiRu) + 

+ (n-l)(n-2) i9ik9lj - 9u9kj)R - 

The crucial fact related to our considerations is that the metric g is conformally flat if and 
only if in the dimension n = 3 the Cotton- York tensor vanishes and in the dimension n = 4 
the Weyl tensor vanishes, see p. 92] or [HJ EH 12] 
Above, Rijki is the Riemannian curvature tensor, 

d d 

jkl ~ dx k jl dx l jk & pk * k pl ' W ~ 9 ijkl 
where T l - k are Christoffel symbols, 

V i 1 n pj ( d 9j P 9g kp dg jk 

jk 2 y \ dx k t- dxi dxP ), 
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Rij is the Ricci curvature tensor, i?^ = R^ k , and R is the scalar curvature R = g^Rij. 
Finally, above V k is the covariant derivative that is defined for a (0,2)-tensor A a and a 
(0,l)-tensor B\ by 

V*4i = - Yl x A vl - r p ki A lp , V*5, = A S| - T p kl B p . 
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